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Abstract. We study spectral variability of 11 ultraluminous X-ray sources (ULX) using archived 
XMM-Newton and Chandra observations. We use three models to describe the observed spectra; a 
power-law, a multi-colour disk (MCD) and a combination of these two models. We find that out of 
the 1 1 ULXs in our sample, 7 ULXs show a correlation between the luminosity and the photon index 
r (hereafter L — T correlation). Furthermore, out of the 7 ULXs that have the L — T correlation, 4 
ULXs also show spectral pivoting in the observed energy band. We also find that two ULXs show 
an L — r anti-correlation. The spectra of 4 ULXs in the sample can be adequately fitted with a 
MCD model. We compare these sources to known black hole binaries (BHB) and find that they 
follow similar paths in their luminosity-temperature (hereafter L — T) diagrams. Finally we show 
that the 'soft excess' reported for many of these ULXs at ~ 0.2 keV seem to follow a trend L cx 
when modeled with a power-law plus a 'cool' MCD model. This is contrary to the expected i cx 
relation that is expected from theory and what is seen for many accreting BHBs. 
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1. INTRODUCTION 

Ultraluminous X-ray sources (ULX) are point like, non-nuclear (not an active galactic 
nucleus, AGN) objects with X-ray luminosities exceeding the Eddington limit for stellar 
mass black holes (StMBH) (Lx ^2 • 10^^ erg s"^). Many ULXs have been identified as 
accreting black holes because of their short- and long-term variability ili and references 
therein]. More than 150 ULXs have been discovered in nearby galaxies 10, [1] and all 
types of galaxies have ULXs L4J- Although they were discovered with the Einstein 
mission almost thirty years ago [seeH, for review], the nature of their X-ray emission 
mechanisms is not yet understood. 

By assuming that ULXs are isotropic emitters, the most luminous source M82 X- 
1 has a luminosity of L ~ iC^^ergs^^ L6J. Therefore, if ULXs do not exceed the 
Eddington limit LEdd ~ 1-3 • 1O^^(M/M0)ergs~^, the accreting black hole should be 
an intermediate-mass black hole (IMBH) with M>1OOOM0. One argument used to 
support the idea of IMBHs powering ULXs is the presence of a ~ 0.2 keV 'soft excess' 
in the spectra of several ULXs [7, 8]. It is not clear, however, whether this excess really 
is a signature of an accretion disk or, maybe, an artifact of a complicated absorber or 
emission from some thermal plasma surrounding the source. Other characteristics seen 
in ULXs in favor of IMBH interpretation (and used against the interpretation of beaming 
of radiation in ULXs [see however 9]) are the low-frequency qua si-p eriodic oscillations 
detected from ULXs M82, Holmberg IX and NGC 5408 lfia[lll.[l2ll. 



There are also several arguments against the hypothesis that ULXs are powered by 
accretion onto IMBHs. The main issue is the formation of such massive objects in 'nor- 
mal' stellar evolution. Black holes of masses greater than 20Mq are not expected to form 
even from the most massive stars ITbI] . Furthermore, about 30% of ULXs are associated 
with star formation regions flT] and a small number of them have been associated with 
short lived, luminous optical counterparts [4]. These facts can be naturally explained if 
ULXs are powered by accretion onto StMBHs. 

There are numerous models put forward on how the luminosity of a StMBH exceeds 
its Eddington limit. Truly super-Eddington models include the advective 'slim' disk 
models 111, [3, 13] or accretion disks with strong outflows |i8iS[i9']. Also the photon 
bubble models poL I21I. I22I] can produce the high observed luminosities. Other plausible 



scenarios involving StMBHs include beaming of radiation by a relativistic jet |23i|24] or 
geometrical beaming by an outflowing wind [|25l l26]. Because of the limited observing 
band of the current X-ray observatories capable of detecting ULXs (~ 0.1 — 10 keV), it 
is not surprising that many of the above models have been successfully used to produce 
the observed spectra of some particular ULX. This is mainly because many previous 
studies have been concentrated only on modelling single observations of each ULX [e.g. 
23,0,121]. 



There are some studies where the variability of ULX spectra have also been studied 
[|30l [31I [32I1 but the samples of ULX with multiple observations has been small. The 
study by Feng and Kaaret [30] showed that the two ULXs in NGC 1313 galaxy have 
different types of spectral variability. Their subsequent study [|33i1 also showed that the 
'soft excess' in NGC 1313 ULX-2 didn't follow the expected L cx relation that is 
expected from theory when modeled with a MCD. 

Therefore, in order to study the spectral variability of ULXs in better detail, we have 
selected a sample of 1 1 ULXs that have been observed more that 5 times with XMM- 
Newton and Chandra. In this paper, we present the results of three simple models that 
we used to describe the observed spectra; a power-law, a multi-colour disk (MCD) 
and a combination of these two models. Our aim is to find correlations between the 
model parameters and similarities between the sources in order to answer the following 
questions: 1) Are ULXs powered by accretion onto IMBHs, StMBHs or are ULXs a 
heterogeneous class of objects? 2) What is the nature of the 'soft excess'? 3) What 
processes produce the observed spectra? 



2. OBSERVATIONS AND DATA REDUCTION 

We searched the data archives for high S /N XMM-Newton and Chandra observations of 
ULXs that had been observed more than 5 times. We also required that the observations 
were not effected by photon pile-up. There were 1 1 sources that met this requirement; 
NGC 253 X-2, NGC 253 X-4, NGC 253 X-9, IC 342 X-6, NGC 1313 ULX-1, NGC 
1313 ULX-2, Holmberg II ULX-1, Holmberg IX ULX-1, NGC 5204 ULX-1, NGC 
5408 ULX-1 and NGC 6946 X-6. In the cases where sources have multiple names, we 
followed the naming convention of Liu and Bregman [3]. 

XMM-Newton data reduction was done using the SAS version 7.1.0. In order to pro- 
duce a homogeneously calibrated data set, we reprocessed all the observations using the 



observation data files (ODF) with the latest calibration files as of March 2008. The reduc- 
tion was done according to the User's guide to the XMM-Newton Science Analysis Sys- 
temQ Specifically, we used the event selections "FLAG==0 && PATTERN<=4" and 
"#XMMEA_EM && PATTERN<=12" for the EPIC-pn and EPIC-mos instruments, 
respectively. Chandra observations were reduced using the CIAO 4.0.1 software. We 
reprocessed all the observations with calibration files from CALDB 3.4.3 following the 
Chandra science threads H 

We used circular source extraction regions of varying radii (between 20" to 40" for 
XMM-Newton and 4" to 12" for Chandra) that covered the full PSF of the sources. Close- 
by source free region^ were used as a background. The resulting background subtracted 
spectral data were binned to require at least 20 counts per bin to ensure adequate statistic 
for the XSPEC spectral fitting. To ensure the best possible spectral quality, the time 
periods of high background were eliminated from the data for both XMM-Newton and 
Chandra observations. 

3. SPECTRAL ANALYSIS AND RESULTS 

We used three models to describe the observed spectra; a power-law, a MCD (diskbb 
model in XSPEC) and a combination of these two models. The effect of the interstellar 
medium was taken into account with the vv'ABS model in XSPEC. The hydrogen column 
density was let as a free parameter but was required to have at least the galactic value 

a. 

All the parameter errors in the following figures are given at the 90% confidence level. 
The errors of the intrinsic luminosities were estimated by repeatedly setting one of the 
free parameters into their minimum (or maximum) values and then fitting the rest of 
the free parameters. We then set the absorption column to zero, computed the intrinsic 
luminosity, and selected the minimum and maximum luminosities as the estimate for the 
error. In the case of the power-law, the luminosities were calculated for the 0.4 — 10 keV 
band. The MCD component luminosities are bolometric. 

3.1. Power-law type ULXs and the L — T correlation 

In Fig. [H we show the results of fitting with the power-law model. There is a corre- 
lation between the luminosity and the photon index T (hereafter L — T correlation) for 
most of the ULXs in our sample. This was also observed for NGC 1313 ULX-1 by Feng 
and Kaaret \M- The ULXs with the L - T correlation are NGC 253 X-4, IC 342 X-6, 
NGC 1313 ULX-1, Holmberg II ULX-1, Holmberg IX ULX-1, NGC 5204 ULX-1 and 



Guide is available at 

http://xmm.esac.esa.int/extemal/xmm_user_support/documentation/sas_usg/USG.pdl 

Threads are available at http://cxc.harvard.edu/ciao/threads/index.html 
^ In the case of XMM-Newton, the background region was chosen from the same EPIC-pn CCD chip as 
the source. 
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FIGURE 1. The intrinsic 0.4 — 10 keV band luminosity versus the photon index. Various observations 
of a given objects showing spectral pivoting in the observed band are connected by a theoretical curve ([TJ. 
Some sources show correlation, but no spectral pivoting and some other even show anti-correlation. NGC 
1313 ULX-1, Holmberg II ULX-1, Holmberg IX ULX-1, NGC 5204 ULX-1 show L-T correlation with 
spectral pivoting in the observed band. L — T correlated sources without spectral pivoting are NGC 253 
X-4, IC 342 X-6 and NGC 5408 ULX-1. NGC 253 X-2 and NGC 1313 ULX-2 are L-T anti-correlated 
sources. NGC 253 X-9 shows no L — T correlation and NGC 6946 X-6 is not variable enough to determine 
any correlations. 



NGC 5408 ULX-1. In contrast, two ULXs of the sample, NGC 253 X-2 and NGC 1313 
ULX-2, show an anti-correlation between the luminosity and T (see Fig. [T] and [30, fig. 
2]). We also find that the luminosity of NGC 253 X-9 is clearly not correlated with T 
and that NGC 6946 X-6 is not significantly variable to determine a possible correlation. 

In general, all the sources that have luminosities permanently above ~ 3 ■ lO^^ergs"^ 
have power-law type spectra. In contrast, the three sources that have luminosities perma- 
nently below ~ 10'^^ erg s"^ have non-power-law (or thermal) type spectra. An exception 
to this 'rule' is NGC 1313 ULX-2 which has non-power-law type spectra at high lumi- 
nosities and power-law type spectra at lower luminosities very much like the very high 
state (VHS) spectra of XTE J 1550-564 L35J- 

Out of the 7 ULXs with the L - T correlation, 4 ULXs (NGC 1313 ULX-1, Holmberg 
II ULX-1, Holmberg IX ULX-1 and NGC 5204 ULX-1) show also pivoting in their 
spectral variability with pivot energies ~ 9,3.5,10,5.5 keV, respectively. Similar 
pivoting has been observed e.g. in Cyg X-1 [36], and several AGNs [see e.g. Isvl table 
1, and references therein]. In Fig. [T] we also plot the dependence between the luminosity 
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FIGURE 2. Left: All observed EFe spectra of NGC 5204 ULX-1 fitted with an absorbed power-law 
model. The curves show the model spectra. Spectral pivoting occurs at ~ 5.5 keV. Right: EFe spectra 
of NGC 253 X-2 fitted with an absorbed power-law plus MCD model. EPIC pn, mosl and mos2 data 
are shown in blue, red and green, respectively. The absorbed MCD component is shown with a pink dot- 
dashed line and the dashed purple line shows the MCD component corrected for absorption. A simple 
absorbed MCD fit with Tin ^ 1.5 keV improves the fit from x^d ^ 1-17 to Xrcd 1-08. Clearly for this 
observation the highly luminous 'cool' ^0.1 keV MCD component results from the improper modeling 
of the underlying power-law component, which leads to almost an order of magnitude increase in the 
modeled absorption column density. 



and the photon index using the expression from [SVi eq. 7] 



2-r 



Le,-e, = AtiD'Fe.-e, = CEl ^ , , (1) 



where D is the distance, Ei and E2 are the energy band values, Fe^-E2 is the flux in 
the selected energy band and C is a constant. Eq. ([T]) give a good representation of 
the observed L — V relation for all the pivoting ULXs. However, NGC 1313 ULX-1, 
Holmberg II ULX-1 and NGC 5204 ULX-1 also show deviations from this when the 
spectrum gets very soft, i.e. when r>3 keV. 

The process producing the observed power-law type spectrum is most likely Compton 
upscattering of soft disk photons in a hot plasma surrounding the disk. The observed 
L — V correlation and the pivoting in the sources can also be explained by variability of 
the luminosity of the soft disk photons. A larger disk luminosity leads to a softer and 
more luminous X-ray spectrum [37]. An example of such variability is for a constant 
hot-plasma luminosity is given in fig. 3 of Zdziarski and Grandi PSQ for 3C 120. The 
similarity of this source to NGC 5204 ULX-1 (and the other 3 ULXs with pivoting) 
is remarkable (see Fig. (2]). Unfortunately the limited observing band of XMM-Newton 
and Chandra do not allow us to see the high energy cutoff of the spectrum which could 
confirm the assumption of the constant hot-plasma luminosity. The deviations from the 
L — r correlation by the three pivoting ULXs can be also understood by a change in the 
mass accretion rate [see [3^, fig 14c]. 



3.2. Non-power-law type ULXs 



Out of the selected 11 ULXs, four ULXs can be adequately fitted with an absorbed 
MCD model. These sources include the two ULXs that show L — T anti-correlation 
(NGC 253 X-2 and NGC 1313 X-2) and the two low luminosity ULXs (NGC 253 
X-4 and NGC 253 X-9) of the sample. The resulting fits are shown in a luminosity- 
temperature (hereafter L — T) diagram in Fig. |3] together with the observations on 
nearby known black hole binaries adopted from Gierlihski and Done [|39i] . Clearly, these 
four ULXs follow the behaviour seen in BHBs. The properties of the sources can be 
summarized as follows. 

(i) NGC 253 X-2 closely follows the L cx T"^ relation. Its track on the diagram is a high 
luminosity-high temperature extension to XTE J1550-564 and LMC X-3. Its relatively 
high luminosity compared to BHB and modeled hot inner disk temperatures could be 
explained by assuming that this source is a maximally rotating Kerr black hole with 
M ~ 7OM0 and inclination of ~ 70 degrees [40] . 

(ii) NGC 253 X-4 on the other hand behaves like LMC X-1 which clearly does not 
follow the expected L oc relation. As discussed in Poutanen et al. [[l9l] . this can be 
understood if the sources are very low mass black holes or even neutron stars accreting 
at slightly above their Eddington accretion rates. In this case the observed temperature 
corresponds to the emission from the spherization radius or even the distant photosphere 
if observed edge-on. 

(iii) NGC 253 X-9 first seems to follow the L cxT'^ track at lower luminosities, but 
then at L > lO'^^ergs^^ deviates towards lower temperatures similarly to LMC X-1 and 
NGC 253 X-4. This type of behaviour is also seen in the spectra of XTE J 1550-564 
when it is in the VHS instead of the high/soft state [35*] . 

(iv) NGC 1313 X-2 shows a slight deviation from the L oc track much like 
XTE J 1650-500. As the luminosity of this source greatly exceeds the Eddington limit 
for a 10 solar mass black hole, the simple diskbb model is not anymore a good 
approximation of the spectrum. Indeed the diskbb model failed to properly model the 
high-energy curvature of the spectrum, predicting slightly too steep cutoff. Because at 
lower luminosities the observed spectra of this source is very much like VHS spectra 
of XTE J1550-564 [|35|l . the higher luminosity non-power-law type spectra could be 
understood as a 'superluminous' state. However, a use of a more elaborate model than 
the DISKBB should be used to investigate this state further. 



3.3. The 'soft excess' and the nature of ULXs 

We show a luminosity-temperature relation for the 'soft excess' in the right panel of 
Fig. [3l Here the model is a combination of a power-law and the diskbb model. We 
find that the intrinsic luminosity of the MCD component roughly scales as L oc T;""^. 
However, this result should be taken with some caution. The highly luminous low 
temperature excesses (L> 10"^° erg s^^. Tin ~ 0.1 keV) do not seem to be real features 
but rather artifacts of improper modeling of the harder emission with a power-law. This 
results in a rise in the modeled absorption column which then causes the increase in 
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FIGURE 3. Left: The luminosity-temperature diagram for the 'non-power-law type' ULXs. All the 
sources follow the tracks of known BHBs. Right: The luminosity-temperature diagram of the 'soft excess' . 
All the sources that can be modeled with the power-law plus cool MCD model follow roughly a i oc Tj^^ 
scaling. Note however that the coolest and highest luminosity 'soft excesses' are artifacts of an improper 
modeling of the hard emission (see Fig.|2]i. 



the estimated intrinsic luminosity of the 'soft excess' (see Fig. O. Never the less, the 
excesses below L < 10^° erg s"^ , Tin ^ 0. 1 keV seem to be real features, but the observed 
trend is not consistent with the interpretation of an accretion disk around an IMBH. 

If the observed soft excess does not originate from a disk around an IMBH, what 
is then its origin? One possible answer could be the 'super-Eddington outflow' model 
by Poutanen et al. [I 19| ] which predicts that the emission from the spherization radius 
becomes cooler with increasing luminosity. However, the observed L oc Tj~^ relation 
is too steep for the predicted relation of this model. The ULXs that have this 'soft 
excess' in our sample have a power-law type spectra. This spectrum most likely arises 
from Comptonization of the disk photons in a surrounding corona. Therefore, a natural 
explanation to both the observed cool accretion disks and power-law type spectra could 
be the 'coupled disk-corona' model [41, i^]. In this model the corona drains power from 
the underlying disk which then effectively cools and becomes less luminous. However, 
there is one important question that this model fails to answer. Assuming that the ULXs 
are at ~ Eddington luminosity sources with masses ~ IOOM0, why would the sources 
be permanently in this 'coronal dominant' state? Or in other words, why are they never 
in a disk dominated (high/soft) state with inner disk temperatures of ~ 0.5 keV? 

Therefore, it is very plausible that many of the above mechanisms operate simulta- 
neously and therefore the ULX variability results from an interplay between outflows, 
advection, beaming (geometrical or relativistic) which is then affected by Comptoniza- 
tion. One possible scenario explaining the nature of the observed spectral variability 
could be a 'hybrid' model combining the 'outflow' model and the 'coupled disk-corona' 



model. At high inclinations the central Comptonized part of the disk would be obscured 
by the outflow. This would result in a non-power-law looking spectra that we see for 
NGC 1313 ULX-2 at high luminosities and for all the NGC 253 ULXs. On the other 
hand, at low inclination angles the central hot corona would be visible to the observer 
at all times and the outflowing wind would cause the Comptonized emission to be geo- 
metrically beamed, thus resulting in a high ~ 10''° ergs~^ apparent luminosity. However, 
it could be evenly plausible that we see power-law and non-power-law type ULXs and 
different kind of variability simply because ULXs are a heterogeneous class of objects. 

4. SUMMARY AND CONCLUSIONS 

We have studied spectral variability of 1 1 ULXs and can draw several conclusions. 

• The ULXs of this sample can be roughly separated into two types based on their 
spectra: power-law- and non-power-law type ULXs. All the power-law type ULXs 
are constantly more luminous than ~ 3 • lO^^ergs"^, i.e. they show permanently 
super-Eddington luminosities for a ~ IOMq black hole. 

• Six out of the seven power-law type ULXs show a L — T correlation, i.e. the 
observed spectrum softens as the luminosity increases. Furthermore, out of the 
7 ULXs with the L — r correlation, 4 ULXs show also pivoting in their spectral 
variability with < 10 keV. Also one non-power-law type ULX (NGC 253 X-4) 
show the L — r correlation. Interestingly, this source is the only one of the non- 
power-law type ULXs that does not follow the L oc T"' track when modeled with a 
MCD. 

• The non-power-law type ULXs can be adequately modeled with a MCD. These 
sources also follow the behaviour of known BHBs in the L — T diagrams. 

• NGC 1313 ULX-2 is a very interesting source showing a VHS type spectra (sim- 
ilar to XTE J 1550-564) at lower luminosities and a non-power-law (thermal) like 
spectra at high luminosities. It is the only source in this sample showing transitions 
between a spectral state that is seen to known BHB and a 'hyperluminous' spectral 
state. 

• The 'soft excess' seems to follow a L oc T^^ relation. Furthermore, the highly 
luminous low-temperature excesses (L> 10''°crgs~^, Tin ~ 0.1 keV) do not seem 
to be real features but rather artifacts of improper modeling of the harder emission. 
This L oc T relation does not support the idea that the observed soft excess arises 
from the disk emission around IMBHs. 

The main question remained to be answered is the origin of these two ULX types or 
classes. It is very plausible that the two observed types arise simply because there are two 
different emission mechanisms producing high luminosities in ULXs. If the ULXs have 
outflows which cause beaming of radiation from the inner parts of the accretion disk, 
the two different kinds of variability can also result from different inclination angles 
between the two types. 
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